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ABSTRACT 

We present the results of a search for molecular gas emission via the CO line 
in the far outer disk of the nearby spiral, NGC 6946. The positions targeted were 
chosen to lie on or near previously- identified outer disk HII regions. Molecular 
gas was clearly detected out to 1.3 R25, with a further tentative detection at 
1.4 R25. The CO detections show excellent agreement with the HI velocities and 
imply beam-averaged column densities of 0.3 — 9 x 10^° cm~^ and molecular gas 
masses of (2-70) x 10^ Mq per 21" beam (560pc). We find evidence for an abrupt 
decrease in the molecular fraction at the edge of the optical disk, similar to 
that seen previously in the azimuthally-averaged areal star formation rate. Our 
observations provide new constraints on the factors that determine the presence 
and detectability of molecular gas in the outskirts of galaxies, and suggest that 
neither the HI column, the metallicity or the local heating rate alone plays a 
dominant role. 

Subject headings: galaxies: individual (NGC 6946) — galaxies: individual (NGC 1058) 
— galaxies: spiral — galaxies: ISM - ISM: molecules - radio lines: ISM 

1. Introduction 

It has become increasingly apparent in recent years that star formation is not just con- 
fined to the optically-bright parts of galaxies. Signposts of low-level star formation, such as 
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HII regions and UV-bright clusters, have been disco vered in the far outer r eg ions of galac 



Ferguson et aLl (Il998a[): iThilker et al. 



Ryan-Weber et al. 



tic dis ks, well beyond the classical R25 radius (e.g. 
( 2005 )) , and in the t idal H I filaments of some interacting systems (e.g. 
( 120041 ) : iBraine et al.l (1200 ll )). The existence of star formation in these unusual environments, 
characterised by low metallicity, low HI column density and low interstellar radiation field, 
provides important constraints on star formation in galaxies, as well as the necessary condi- 
tions for sustaining a multi-phase interstellar medium (ISM). 

Much theoretical effort has been devoted to understanding the formation of a cold 
molecular phase in the ISM, a prerequisite for star formation. Thermal instabilities are 
generally believed to be the primary mechanism for converting cool atomic gas into cold 
molecular gas once a mmimum pressure or local threshold surface density is exceeded (e.g. 
Elme green fc Parravand (jl994j ) ; ISchayd (120041 ) ; iBlitz fc RosolowskyI (120061 ) ) . Additional mech- 
anisms for forming unstable cold clouds includ e shocks frorn spiral density wav e s and swing- 
amplifier and magneto- Jeans instabiities (e.g. iDobbs et al.l (120061 ): iKim et al.l (120021 )). The 
threshold column densities predicted for molecular cloud formation are on aver age a few 
times s maller than those in f erred for massive star formation in galaxies (e.g. 
(1l987l ): iMartin k KennicuttI J200l[ )l 



Skillman 



Direct detection of molecular clouds in extreme environments provides a means to test 
ideas about cloud formation and survival. Carbon monoxide (CO) is the most abundant 
heteronuclear molecule (i.e. with permitted rotational transitions) and is excited by collisions 
at densities typical of molecular clouds, m aking CO the standa rd tracer of molecular gas. 
Concerns have previously been raised (e.g. iPfenniger et al.lll994l ) as to whether the physical 
conditions in the outskirts of disks are sufficient to excite CO to detectable temperatures 
and whether the metallicities in these regions are sufficient to form enough CO for detection. 
With the discovery of CO emission out to 1.5 R25 in the disk of NGC 4 414 (IBraine fc Herpin 



2OO4J) and in high column density tidal HI filaments in the M81 Group (IBrouillet et al. 



11992 



Walter et al.ll2006l ). these concerns have been proven somewhat unfounded. 



In this Letter, we present the results of a search for molecular gas in the far outer disk of 
the nearby spiral N GC 6946. While the stu dy of molecular gas in NGC 4414 focused on high 
column HI regions (IBraine &: Herpiru 120041 ). we have chosen to search for CO emission in the 
vicinity of outer disk star-forming regions. NGC 6946 was previously identi fied through Ha 
observ ations to exhibit very extended massive star formation, out to ~ 2R25 ( iFerguson et al. 



1998al ). and thus presents an excellent opportunity to attempt to measure and quantify 
molecular gas at extreme radii. We also report on the null detection of molecular gas at two 
positions in the far outer disk of another spiral with extended star formation, NGC 1058. 
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2. Observations 

Our search was conducted using tlie IRAM 30m telescope at Pico Veleta, Spain during 
July 7-13 1999 in very good summer weather. The ^^CO J= 1 — and J= 2 — 1 transitions 
at 115.271 and 230.538 GHz were observed simultaneously in two polarizations (H and V) 
with the newly installed "AB" SIS receivers. The filter bank back end was configured into 
four units of 256 x IMHz channels, one for each polarization and each frequency, yielding a 
650 km s"'^ bandpass and a channel separation of ~ 2.6 km s~^ in the CO(l-O) transition, 
respectively 330 and 1.3 km s"^ in the CO (2-1) transition. Typical system temperatures 
(T;) for the observations of NGC 6946 and NGC 1058 were 250 K at 115 GHz and 400 K at 
230 GHz. All spectra are presented on the main beam temperature scale assuming forward 
and main beam efficiencies of 0.9 and 0.73 at 115 GHz and 0.84 and 0.51 at 230GHz. 

Single pointings were obtained on or near Hll regions lying at radii in the range 0.7- 
1.6 R25 in NGC 6946 (see Figure [T]), and 1.0-1.9 R25 in NGC 1058. The IRAM 30m beam 
size is 21" at 115 GHz and 11" at 230 GHz , corresponding to 560 and 280 pc respectively 
at the distance of NGC 6946 (5.5 Mpc, see Regan et al. ( 200^ ) and 1 kpc and 0.5 kpc at 



the distance of NGC 1058 (10 Mpc, see iFerguson et al.l (jl998al )). As can be seen in Figure 



[H the CO beam size was usually larger than the size of the target HII regions. 

After eliminating obviously bad spectra and interpolating over occasional bad chan- 
nels, spectra were summed for each transition and position. Despite observing in position- 
switching mode with a reference position typically ~ 3 R25 from the galaxy, only a zero-order 
baseline subtraction was required. 



3. Results 

The results of our CO search are presented in Table 1. In NGC 6946, CO is clearly 
detected at 7 out of 10 positions, ranging from 0.7-1.3 R25, with a further tentative detection 
at I.4R25 (P3). In NGC 1058, no detections were made. 

Figure [2] shows the CO(l-O) (and C0(2-l) when detected) spectra of the eight positions 
in NGC 6946 with molecular gas emission. CO (2-1) emission is detected at most positions 
with CO(l-O) emission which suggests that the gas is not highly subthermally excited, as it 
might be if it were extremely diffuse. HI spectra extracted at the same position s and with 



the sa me synthesized beam size are overlaid, taken from the 21cm data cube of iBoomsma 



(120071 ). The HI and CO velocities show excellent agreement, with the maximum offset being 
less than 3 km s~^. The CO line widths, as determined from gaussian fits, are in the range 
8-30 km s~^ and are greater than those expected from individual Giant Molecular Clouds 
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(GMCs), even at the outermost positions. This imphes that an ensemble of clouds is observed 



within the 560 pc CO(l-O) beam in each case, similar to the findings of iBraine fc Herpin 



fl2004h for the outer disk of NGC 4414. 



The spectra at the two outermost positions warrant further discussion. The systemic 
velocity of NGC 6946 is very close to the velocity of the HI and CO in the Milky Way, thus 
making comparisons at some positions more difficult than they would be in higher velocity 
systems. The detection at P8 (1.3 R25) could be called into question since the emission line 
is separated by only a single velocity channel from the narrow Galactic line at this position. 
If the line emission were Galactic in origin then we would expect it to have a significant 
spatial extent on the sky and to be detected at positions P5 and P4. For a CMC lying 
500 pc along the line of sight (i.e. 100 pc above the Galactic plane), P5 and P4 correspond 
to linear distances of 0.12 and 0.16 pc from P8. No CO emission at the velocity detected 
at P8 is seen at these positions. Fig. 2 shows that weak HI is detected at the local velocity 
(~ — 15 km s"^) by the WSRT in all positions but there is likely some spatial filtering out 
of the local HI signal. At PS the HI is much stronger at the velocity of the CO we attribute 
to NGC 6946, further strengthening our argument that the P8 detection is real and that the 
line belongs to NGC 6946. We thus consider P8 to be a genuine detection of CO emission 
from NGC 6946. 

The line emission at 150 km s~^ towards P3 (1.4 R25) is statistically significant (~ 5a) 
but another "line" of similar strength is seen at 270 km s~^. The precise agreement between 
the CO and HI line velocities at this position certainly supports the reality of the CO feature 
but the presence of the other peak casts some doubt. As a result, we consider P3 as a probable 
but still tentative detection. 

Figure [3] shows a view of the spiral arm segment containing positions P2, P7, and P9- 
P12 in the Hct line and in 8/im emission, usually attributed to PAHs, with HI contours 
overlaid. The PAH emission was observed as part of the SINGS project with Spitzer. While 
the locations of the HI, Ha, and PAH 8/im emission are closely correlated, the intensities of 
the emission at different wavelengths do not correlate. While Pll and P2 have similar Ha 
intensities (see also Table 1), the 8/im and CO emission is much weaker further out. The 
N(H2)/N(HI) ratio also decreases, by more than an order of magnitude, from 0.7 to 1.4 R25 
and is largely driven by a sharp decrease in the CO intensity as one crosses R25 (P12 to P7). 

Table 1 also lists the implied H2 column density at each position, derive d using a "stan- 



dard" conversion factor of N(H2)/Ico = 2 x 10^° H2 cm ^ per K km s ^ (jPickman et al. 



19861 ). and the molecular gas mass within the beam, ir icluding a correcti on for Helium. The 



masses inferred are in the range ~ 2 - 70 X lO^M 0. iDigel et all fll994f ) studied molecular 



clouds in the far outer disk of the Milky Way (presumably beyond R25 assuming a standard 
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rotation curve) and found clouds of much lower mass, in the range 2 — 40 x lO^M q. If such 
clouds also exist in the outskirts of NGC 6946, they would individually lie well below our 
detection limit of ~ lO^M fr,. 



Discussion 



We have presented evidence for the existence of molecular gas in the far outskirts of 
NGC 6946 and noted a rather abrupt decrease in the molecular gas fraction, as well as the 
PAH emission, as one crosses R25 along an individual spiral arm. An important question is 
whether this behavior could be explained by metallicity variations in the unde rlying gas disk. 



Indee d there is strong evidence for a radially increasing N(H2)/Ico v alue (e.g. ISodroski et al. 



19951 ) as well as much higher values in low-metallicity systems (e.g. iRubio et al.lll99ll ). The 
metal licity in the outer disk of NGC 6946 varies by a factor ~ 5 from log(0/H)= —3.2 to 
—3.9 (IFerguson et al.lll998bl ) while the Ico/N(HI) ratio decreases by a factor of ~ 30. Focus- 
ing only on the positions which bridge R25 (P12 to P2), the metallicity remains essentially 
constant while the Ico/N(HI) ratio decreases by a factor of ~ 6. In the absence of a very 
strong rad iation fi e ld, th e N(H2)/Ico factor should not vary more than linearly with metal- 
licity and IWilsonI (119951 ) estimates that N(H2)/Ico oc from a study of molecular 
clouds in Local Group galaxies. It is thus reasonable to expect that the N(H2)/Ico ratio 
might be a factor 2-3 higher in the most distant points than at P9, where the "standard" 
value is probably appropriate. Such a correction is not sufficient to explain the strong de- 
crease in Ico/N(HI) at R25 and we conclude that there is a genuine drop in the molecular 
gas content at the edge of the optical disk as defined by R25. 

Although our pointings were selected to lie on or near outer disk HII regions, the Ha 
luminosity contained within each beam shows significant variation and co rresponds t o loca l 
star formation rates (SFRs) ranging from 10~^ — 10~^ Mq yr~^, using the iKennicuttI (119981 ) 
calibration (SFR=L(Ha;/1.26 x 10^^ erg s~^). The Ha emission in some regions is predom- 
inantly diffuse (e.g. P12 and P7) and may more refiect ionizat ion by photons produced 
elsewhere rather than in situ star formation (IFerguson et al.lll996l ). The star formation effi- 
ciency (SFE), defined as SFR/M(H2), ranges from ~ 10 — 0.07 Gyr~^; eliminating the two 
outlier points (P12 and P2), the range reduces to 2 — 0.5 Gyr~^. The mere presence of 
a luminous HII region does not guarantee the detection of molecular gas. Position PI in 
NGC 1058 is the second most luminous star-forming region targeted in our study, and is 
our second most sensitive integration, yet CO was not detected there. Intriguingly, there 
is a potential correlation between the detection of CO (2-1) emission and the presence of a 
luminous HII region in that the three non-detections of CO (2-1) are among the four lowest 
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Ha fluxes. 



With the exception of P4, the HI column densities at our observed positions within 
NGC 6946 are quite high, close to or greater than lO^^cm"^ on the scale of several hundred 
parsecs. The HI column density may be an important factor in determining whether molec- 
ular gas will be det ected but a high HI column is clearly not a sufficient condition (see also 
Gardan et allboOTh . 



It may be purely coincidental that N(H2) /N(HI) changes abruptly across R25 in NGC 6946 
While R25 provides a means to characterise the optical extents of galaxies, it is not expected 
to relate to any underlying physical properties. In NGC 6946, neither the metallicity, the HI 
surface density nor the surface brightness of the stellar disk exhibit unusual features at this 
location. Interestingly, the azimuthally-averaged areal star formation rate in NGC 6946 does 
show a sharp decline near the edge of the optical disk (~ 0.8 R25), although widespread low- 



level star formation is observed much further out (IFerguson et al.lll998al : lMartin fc Kennicutt 



Elmeereen & Parravano 


1994 


Jschay. 


2004) 



although the physics that dr i ves this behavior is still open to debate. Pressure based 



( Elnaegreen &: Parravandll994l : iBlitz fc Rosolowskyll2006l ) and column density based (ISchaye 
20041 ) models are very difficult to distinguish once the HI surface density exceeds that of the 
stars. Larger-scale mapping of the outer disk of NGC 6946 and sensitive observations of 
the outskirts of additional galaxies are required in order to more thoroughly understand the 
molecular content and distribution at large galactocentric radii. 



AMNF acknowledges a Marie Curie Excellence Grant from the European Commission 
under contract MCEXT-CT-2005-025869 and a Visiting Professorship at the Observatoire 
de Bordeaux where this work was completed. We thank Reuse Boomsma and Tom Oosterloo 
for kindly sharing their HI data on NGC 6946 prior to publication. 
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Table 1. Molecular Gas in the Extended Disk of NGC 6946 



Posn. 


R.A. 


Doc. 


R/R25 


ICO{1-0) 


rms 


AV 


IcO(2-l) 


rms 


(O/H) 


N(H2) 


M(H2) 


N(HI) 


L(Ha) 




(J2000.0) 


(J2000.0) 




K km s-l 


mK 


km s~ 


1 K km s-l 


mK 


dex 


1020cm-2 


lO^MO 


1020cm-2 


lO^^erg s-l 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 




NGC 6946 


P9 


20 34 34.9 


60 11 38.4 


0.7 


4.5±.l 


17 


19 


6.3±.3 


39 


-3.2 


9.0 


70 


22.1 


142.5 


PIO 


20 34 25.9 


60 11 31.3 


0.8 


2.3±.l 


16 


22 


2.2±.2 


30 


-3.3 


4.6 


36 


17.1 


64.3 


Pll 


20 34 25.9 


60 11 51.3 


0.8 


2.2±.l 


15 


24 


3.8±.3 


32 


-3.3 


4.5 


34 


15.2 


28.9 


P12 


20 34 25.9 


60 12 11.3 


0.9 


l.Oit.l 


10 


30 


< 0.5 


23 


-3.4 


2.0 


15 


10.1 


1.4 


P7 


20 34 26.6 


60 12 47.9 


1.0 


0.21±.03 


4.3 


18 


0.22±.05 


7 


-3.4 


0.4 


3.3 


9.2 


2.2 


P2 


20 34 29.0 


60 13 50.5 


1.1 


0.14±.02 


3.0 


14 


0.14±.03 


5 


-3.5 


0.3 


2.2 


8.3 


28.8 


P8 


20 34 54.3 


60 15 40.5 


1.3 


0.15±.02 


5.0 


8 


< 0.3 


13 


-3.6 


0.3 


2.3 


8.9 


4.2 


P3 


20 33 48.4 


60 08 32.6 


1.4 


0.11±.02 


3.7 


12 


< 0.13 


8 


-3.9 


0.2 


1.7 


13.6 


4.7 


P5 


20 34 51.4 


60 16 24.9 


1.5 


< .13 


6.1 




< 0.5 


24 


-3.8 


< .3 


< 1.9 


10.2 


< 0.7 


P4 


20 34 51.4 


60 16 44.9 


1.6 


< 0.14 


6.7 




< 0.4 


20 


-3.8 


< .3 


< 2.0 


5.0 


< 0.1 


NGC 1058 


PI 


02 43 29.3 


37 18 57.6 


1.0 


< 0.07 


3.5 








-3.3 


< 0.14 




3.7 


85.5 


P2 


02 43 29.4 


37 17 40.0 


1.9 


< 0.11 


5.4 








-3.8 


< 0.21 




4.2 


11.7 



Note. — Units of right ascension are hours, minutes and seconds and units of declination are degrees, arcminutes and arcseconds. Col.(l): position ID. 
Cols. (2) and (3): Right ascension and declina tion of the pointing. Col. (4) the dcpro jcctod radius of the pointing in terms of R25, calculated for NGC 6946 
using R25 = 5.75' jde Vaucouleur s et al.| [ l99ll ). and i = 34.0° and P.A.= 69° tConsidc rc fc Athana ssoula 198^) and for the almost face-on NGC 1058 using 
R25 = 1.51' llde Vaucouleurs et al.lll99ll) and i = 0° and P.A.= 0°. Cols.(5) and (8): the CO(l-O) and CO(2-l) intensities on the main beam temperature scale 
respectively. Upper limits are 3(t and calculated assuming a 7 channel (18.2 km s"!) line width. Cols.(6) and (9): the rms noise, measured per 2.6 km s 
channe l. Col. (7): the velocity width as determined from fitting gaussians to the CO(l-O) spectra. Col. (10) The oxygen abundance, from l lFerguson et al.l 
lll998bl) and Ferguson, unpublished. Note that the metallicities for NGC 6946-P3 and NGC 1058-Pl and P2 are directly measured, whereas those for other 
positions are inferred from interpolation of the abu ndance gradient. Cols. (11) and (12) the H2 column density and cloud mass, as described in the text. 
Col. (13) The HI column density at the position from iBoomsm 3 (|2007|) and lOickev et al.l lll990l) . Col.(14) The Ha luminosity within a 21" aperture, corrected 
for Galactic extinction and [Nil] emission. 
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Fig. 1. — An Ha image of NGC 6946. The white and black circles show the IRAM positions 
observed and the size is that of the CO (1-0) beam. The black ellipse shows the R25 contour. 
The box shows the region covered in Figure [3l 
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Fig. 2. — CO spectra of the detections in NGC 6946, along with P3 which is a tentative 
detection. The vertical axis shows the main beam antenna temperature expressed in Kelvins 
and the horizontal axis the heliocentric velocity in km s~^. The CO (1-0) spectra are shown 
as solid black lines, the CO (2-1) as green dashed lines when detected, and the HI spectra 
as red dotted lines. The HI scale is constant with respect to CO in order to show how 
relative HI and CO line strengths vary from position to position. Dividing the CO scale by 
3 gives the HI flux in Jy beam~^. The narrow lines at— 15<^^<— 8 km s^^ are due to 
local emission from high-latitude Galactic clouds along the line of sight (b ~ 12°). Note the 
vertical scale varies between plots. 
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Fig. 3. — (Top) An Ha image of the spiral arm segment containing the P2, P7 and P9-12 
positions, with the CO(l-O) beam superposed. The HI emission is shown as yellow contours 
at the 5, 10, 15, and 20 xlO^° cm~^ column density levels. (Bottom) As above but overlaid 
on the Spitzer IRAC 8 /xm image. 



